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Abstract 

We calculate the 'exact' potential corresponding to a one-dimensional interacting system of two 
electrons with a specific, tailored density. We use one-dimensional density-functional theory with a 
local-density approximation (LDA) on the same system and calculate densities and energies, which 
are compared with the 'exact' ones. The 'interacting-LDA system' [j| corresponding to the LDA 
density is then found and its potential compared with the original one. Finally we calculate and 
compare the spatial entanglement of the electronic systems corresponding to the interacting-LDA 
and original interacting system. 
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I. INTRODUCTION 



The accuracy of the entanglement calculated using approximations is important when 
modelling the suitability of electron systems for use as quantum information devices. Inves- 
tigating this accuracy can also reveal more information about an approximation and perhaps 
ways of improving it. We look at the ground state (GS) of a one-dimensional (ID) contact- 
interacting system of two electrons. The electronic spins are maximally entangled (singlet 
state), so we focus on the spatial entanglement l|. We model the system approximately using 
density-functional theory (DFT)^ with the one-dimensional analogue of the local-density 
approximation (LDA) . DFT gives the GS density, not the many-body wave-function, and 
although in theory all GS properties are expressible as functionals of the GS density, it is not 
known how to calculate the entanglement this way. An alternative route is to find the (con- 
tact) interacting LDA system (i-LDA) [1] that gives the LDA density. We may now access 
the entanglement relatively easily from the i-LDA GS many-body wavefunctionjl]. By the 



Hohenberg-Kohn theorem 



2j the GS density uniquely determines the external potential for 



a given form of interaction, then the i-LDA system must be the unique contact-interacting 



system that corresponds exactly to the LDA density. In J4] we introduced the iterative 
scheme 

<£tVi) = — Url^lkM - n^W] + vUri) (1) 

to find the potential v ex t of the interacting system that reproduced a density arising from 
a single particle equation. This scheme may also be used in a more general way, i.e. to 
find the interacting system that gives an experimental density, or to design a potential to 
give a desired (target) density. The latter application could be important in the future 
for the design of nano-devices. Hence in this contribution we illustrate the usefulness of 
the iterative scheme (Eq. [1]) in finding the numerically 'exact' potential that gives a target 
density in a simple one dimensional two electron system of Hamiltonian H = 5(xi — X2) + 
J2i=i 2 — + v ext(X'i) (atomic units). We also check the scheme by showing that it can 
correctly reproduce a pre-chosen potential from the corresponding density. We then apply 
DFT and the LDA using the 'exact' potential and again use the scheme to show how close 
the i-LDA potential is to the exact and to calculate the approximation to the entanglement 
corresponding to using the LDA. 



To apply the iterative scheme Eq. [TJto the tailored GS density and find the potential for 
the system, we expand the interacting GS using the single particle eigenvectors corresponding 
to f* xt . This basis is changed at every step and the components of the GS are symmetrised. 
This means that the calculations become much longer but tractable for a one dimensional 
system with a contact interaction. By changing the basis at each step we can access a much 
larger set of wave- functions with a small basis size so in theory allowing us to approximately 
reproduce any v-representable density. For fast calculations we use a single particle basis of 
10 thereby giving a symmetric two particle basis of 55. 

We investigate densities on the range [—4 : 4] with 200 mesh points and start with the 
trial potential vl xt (x) = O.OOOlx 2 . This produces a relatively spread out density on this range 
so reduces the risk of the iterative scheme exploding at the start due to the denominator 
being too close to zero. The error is quantified using i Y^i=i \n tar9et (xi) — n tnal (xi)\ where 
iV is the number of mesh points (200 in this case). The potential that gives the density to 
high enough accuracy is then used in the ID DFT scheme to produce the LDA GS density. 



II. ONE-DIMENSIONAL DFT USING LDA 



We use the one- dimensional DFT for contact interacting fermions of Magyar and Burke 
j^l with their local-density approximation. DFT seeks to express the total GS energy as a 
functional of the density n(r), E[n] = T NI [n] + U[n] + E xc [n] + f v ext (x)n(x)dx. This gives 
rise to the single-particle Kohn-Sham equations 

1 d 2 



2dx 2 



+ v ext (x) + v H (r; [n]) + v xc (r; [n]) <pi(x) = e^x 



(2) 



thereby allowing the density n = Y2i |0«| 2 t° be efficiently computed. 

Here vh = W- and v xc = ^ x& . For a contact interaction U = J (n 2 /2)dx so Vjj = n. E xc 
can be written in the form E xc = J e xc (n)n(x)dx where e xc = e x + e c . The results for a 
one dimensional homogeneous electron gas are used to give the local-density approximations 
= -„/4 and e-» = („• + ^ + *. + e), where „ = - 1/2 4, „ = -0.0043 61 43, 
d = 0.252758, e = 0.0174457 and we have used the Pade parametrization of Ref. p(. 

We then take the functional derivative of E xc to obtain v xc = v x +v c . This gives v x (n) = 
—n/2 and 



v c (n) 



l + n 



2an + b 



2n + d 



an 2 + bn n 2 + dn + e 



e c (n). 



(3) 
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We may now use the LDA density with the iterative scheme Eq. [I] to calculate the i-LDA 
potential. This is shifted by a constant to give equality of energy between the i-LDA and 
LDA system. The procedure starting from a density with an unknown potential may be 
summarised as 

^target E 1- HI 'exact' E 1- El LDA E 1- ffl „i-LDA ( A\ 

u u ext u v ext ■ V*) 



III. SYSTEM 1 



We first consider a known potential v 



target 
ext 



X 



0.02x+sin(3x) + 0.3x 2 and the GS density 



n target aT [ s [ n g f rom it. We wish to test the scheme's ability to reproduce the target density 
and potential. 
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FIG. 1: The target density, 'exact' density 
calculated using the 'exact' potential arising 
from the iterative scheme and the LDA den- 
sity obtained using the target potential 



FIG. 2: The target potential and the 'exact' 
and i-LDA potential resulting from the 
iterative scheme. 



We see in Fig. [TJthat the density has been reproduced almost exactly (~ 8.5 x 10 -6 error) 
by using Eq. [T]with n tar9et . We shift the resulting 'exact' potential by a constant so it gives 
the same energy as v^ t 9e and find that the 'exact' potential is almost identical to the known 
potential: they cannot be distinguished by eye (Fig. and differ only at the extremities 
where the density is low. We then use v^ e to calculate the LDA density which is a relative 
good match to the exact density, see Fig. [TJ We now calculate the i-LDA potential which 
gives the LDA density to a high accuracy (~ 8.6 x I0~ 6 error). Fig. [2] shows that the exact 
and i-LDA potentials are very different at the extremities but this is relatively unimportant 
as the density is negligible in this region. An enlarged view of the region where the density 



is non-negligible (inset) shows that the i-LDA is a good match to the exact potential for 
most of the interval. 



IV. SYSTEM 2 



We now apply the iterative scheme to a density for which we do not know the potential. 
We choose a flat-topped density which dies off as a Gaussian for values close to the box 
boundary as the target density. 

We see in Fig. [3] that the target density is reproduced well (~ 5.7 x 10~ 5 error) by the 
'exact' potential shown in Fig. HI This potential is then used to calculate the LDA density 
which is shown in Fig. [3] and is a good match except in the flat region. The i-LDA potential 
gives the LDA density with ~ 1.5 x 10~ 4 error and it is very different from the 'exact' 
potential at the range boundaries. Here though the density is almost zero so the accuracy of 
the potential would not be expected to be high in this region. For areas where the density 
is non-negligible we see that the i-LDA potential is similar to the 'exact'. 
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FIG. 3: The target density, 'exact' density 
using the 'exact' potential arising from the 
iterative scheme and the LDA density using 
this potential 



FIG. 4: The 'exact' and i-LDA potential for 
a density with a flat top. 



V. ENTANGLEMENT 



We may use the i- 
sponding to the LDA 



DA system to calculate the electronic spatial entanglement corre- 
l|. We quantify the spatial entanglement using the Von Neumann 





Exact Energy 


LDA Energy 


Exact Entanglement 


LDA Entanglement 


System 1 


0.832 


0.821 


0.161 


0.140 


System 2 


0.683 


0.692 


0.206 


0.194 



TABLE I: Table showing the energy of the exact and LDA system in Hartrees and the exact 
and LDA spatial entanglement quantified by S the Von Neumann entropy of the reduced density 
matrix. 

entropy of the reduced density matrix S = — Trp re d log 2 p re d which we calculate by express- 
ing p re( j in terms of the basis and diagonalising. Results are shown in Table [H They show 
that the ID LDA energies are very accurate for the densities considered here, corroborating 
the results in Also 'exact' and LDA spatial entanglement are fairly close. This is in- 
teresting as in earlier work[l] approximations to the LDA entanglement for the 3D Hooke's 
atom suggested a much less similar entanglement especially for values of S as large as found 
here. We note though that for a confining potential with strength u = 0.5 the 3D system 
was deemed as low interacting, with entanglement S = 0.0228, and could be modelled well 
by the LDA. However for ID Hooke's atom and u = 0.5 we find S = 0.194 - similar to 
the magnitude of the entanglement found in this paper - suggesting that the densities we 
investigated here also correspond to not too strongly interacting systems. This implies that 
LDA should be accurate for the corresponding entanglement, which supports our findings. 

We have demonstrated that the scheme Eq. [T] can be used to find a potential which 
reproduces a target ground state density of an interacting system. In addition the scheme 
may be used to find the entanglement corresponding to the LDA. 
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